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THEORETICAL LAMINAR CONVECTIVE HEAT TRANSFER 

AND BOUNDARY-LAYER CHARACTERISTICS 

ON CONES AT SPEEDS TO 2k KM/SEC 

By Gary T- Chapman 

Ames Research Center 
Moffett Field, Calif. 

SUMMARY 


Solutions of the equations of Boundary layers without pressure gradient 
for air in chemical and thermodynamic equilibrium have Been obtained for 
free-stream velocities from 3 to 2 h kilometers per second, wall temperatures 
from 1000° K to ^200° K, Boundary-layer pressure levels from 0.1 to 10 atmos- 
pheres, and values of the ratio of the kinetic energy per unit mass of the 
stream at the edge of the Boundary layer to the total enthalpy (i.e., Ue 2 /2He) 
of from 0 to 0-93- From these solutions the heat transfer, skin friction, 
Boundary-layer thickness, displacement thickness, and momentum thickness were 
calculated. Simple correlations were obtained for the heat-transfer and skin- 
friction results. The effects of various flow conditions on these correla- 
tions as well as on the noncorrelated results were assessed. A complete 
tabulation of all the necessary parameters is also included so that additional 
correlations may Be sought . 

Besides the presentation of the solutions a procedure is outlined which 
illustrates how the use of the Crocco equation leads to very good first 
approximations to start the iteration procedure used in obtaining the present 
solutions. The use of the Crocco equation to start the iteration procedure 
results in very rapid convergence to the final solution. 


INTRODUCTION 


The prospect of manned space flight to the planets has stimulated con- 
siderable interest in heat transfer to vehicles returning to earth at speeds 
considerably higher than escape speed. It has been shown (ref. l) that 
radiative heat transfer to blunt shapes can be an important or even the domi- 
nant mode of heating in this speed range. A recent study by Allen, Seiff, and 
Winovich (ref. 2) has shown that the radiative heating and also the total 
heating can be reduced if sharp or slightly blunted cones are used instead of 
bluff bodies . 

Many studies have been made of the convective heat transfer to the 
stagnation region of blunt bodies at high speeds (refs. 3* and 6) and to 

blunt-nosed slender bodies (refs. 7 and 8) through the speed range where 
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dissociation but little or no ionization occurs. However., boundary layers on 
cones have not been treated for speeds where dissociation and ionization 
occur. In reference 2 , the method of reference 7 "was applied to the boundary 
layer on high-speed cones, but this required extrapolation well beyond the 
conditions considered in reference 7. The present study was undertaken to 
provide a more rigorous treatment of the heat-transfer and boundary-layer 
characteristics for cones at still higher speeds, where the effects of both 
dissociation and ionization are important. 

Calculated heat transfer, skin friction, and three boundary-layer thick- 
nesses are presented' for equilibrium flow over sharp cones. These solutions 
are compared with the Crocco equation for real gas; and a procedure is 
developed for using the Crocco equation as a starting point in solving the 
more general boundary-layer equations . 


SYMBOLS 


a ,b ,m 

CDf 

Cp 

Cf 

C P 

H 

h 

K(p e ,T v ) 

k 

l 

Hu 

n 

P 

4 


skin-friction correlation constants 
skin-friction drag coefficient 
average skin-friction coefficient 
local skin-friction coefficient 

total specific heat (including the chemical reaction terms) 

normalized variables defined by equation (8) 

u 2 

total enthalpy, h + — 
static enthalpy 

heat-transfer correlation constant 

total thermal conductivity (including the chemical reaction terms) 
PP 

PwPw 

<i c Pv x 

Nusselt number. r 

(H e - H w )k w 

exponent, 0 for two-dimensional flow and 1 for three-dimensional 
flow 

pressure 

heat-transfer rate 
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I? 


R 

Re 


Rw 


r 

St 

T 

u,v 

x*y 


P 

6 

5* 


^i,cfe,cp 3 ,cp 4 

P 

Q 

0 c 

M* 

T 

o 


gas constant 

Reynolds number, PeUeX 

Pe 

Reynolds number, 

radius of cross section 

Stanton number, 2= 

(He - Hv)peUe 

temperature 

components of velocity in x and y directions 

coordinates along a body meridian and normal to a body 
generator 

pressure gradient parameter (see eq. ( 9 )) 

boundary-layer thickness to — = 0.995 

Ue 

boundary-layer displacement thickness 

convergence criteria (see eqs . ( 15 ) , (l6), ( 17 ) , and (l8) ) 
air density 

boundary- layer momentum thickness 
cone half-angle 
coefficient of viscosity 
shear stress 

transformed coordinates defined by equation (7) 
total Prandtl number, |j, ^ 


Subscripts 


b at model base 

e conditions at edge of boundary layer 

i number of iteration 
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st stagnation conditions 

w wall conditions 

tj , £ differentiation with respect to 

00 free-stream conditions 


ANALYSIS 


The "boundary-layer equations and the initial steps in their solution are 
first stated in general terms in order to illustrate a wider application of 
the method of solution. Subsequently, the simplifications pertinent to flow 
over isothermal cones are made . 


General Equations 


The general equations for a thin laminar boundary layer for a dis- 
sociating and ionizing gas in equilibrium are, in the Coordinate System defined 
by sketch (a), the continuity equation 



d(pur n ) + 5(pvr n ) = Q 


dx 


dy 


the momentum equation 


pu £^ + pv £^=_^P + d[p(du/dy)j 
dx dy dx dy 


and the energy equation 


dh . dh dp /du 

pu r + p v r = u n + ^ r 
ox dy dx \dy 


2 d 


k_ dh 
C P dy 


dy 


( 1 ) 


( 2 ) 


( 3 ) 
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The gas state and transport properties cp, k, and p include the chemical 
reaction terms implicitly (i.e., they are total properties as defined by 
Hansen, ref. 8). The exponent n in equation (l) is 0 for two-dimensional 
flow and 1 for axially symmetric flow. Equations (l) , (2), and ( 3 ) are based 
on the usual thin boundary-layer assumptions . Energy transport by radiation 
has been neglected. The neglect of radiative transport is realistic if the 
ratio of the energy radiated (per unit mass along a stream line) to the local 
static enthalpy is small. The range of applicability of this assumption is 
hard to define because this ratio depends on the free-stream velocity, cone 
angle, and the distance along the stream line, and will for the case of a 
small angle cone have its maximum within the boundary layer. For this case 
the neglect of energy transport by radiation is conservative in that the con- 
vective heat-transfer rates thus calculated are on the high side. The pres- 
sure gradient term has been retained for generality but will be set equal to 
zero for the cone solutions . 

By letting 


H = h + 


and combining equation (2) with equation (3)> we obtain 


dH ^ 3 h 

pu — + pv — 
ox ay 


dfiiM 
9 ay 
dy 


+ — ! !=■ 


1 -i 

cr, 


u= 


dy/J 


Sy 




(5) 


where a is the Prandtl number. 

The solution of the boundary-layer problem with zero mass injection 
requires solving equations (l), (2), and (5) simultaneously for boundary 
conditions 

at y = 0 

u = v = 0 and H = H w | 


and as y -*■ oo 


u -* Uq and H -*■ H e 


(6) 


Equations (l), (2), and (5) may be put into a form more suitable for 
solution by use of the Howarth and Mangier transformations 


Tl 3 r p dy 

7^ J o 


(7) 
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-r- 


5 = J Pw^ u e^ " is 

O 


and the following dimensionless dependent variables 


= f = 1 f n dl »' and g = i§ 


( 8 ) 


Substituting equations (l), (7), and ( 8 ) into equations (2) and (5) yields 


and 


where 


and 


( TJTj) T) + ff 

(i znX + f §n + 


W + 

u e 2 

H e 


2p (^ ' f,)2 ) ' 2?(f i f ^ - f £ f W 
_( 1 - i) - 25(f,g 5 - 

1 = «*/ Pw^ 


( 9 ) 


(10) 


_ d in tie 
d In ^ 


The boundary conditions transform to 
at t] = 0 


f = ^ = 0 g = Sw 

and as 7 ] oo (jjJ 

^Tj "*• 1 g -*• 1 

* 

Gas properties .- The solution of equations ( 9 ) and (10) requires an 
equation of state and equations for the transport properties . As these 
properties are allowed to vary, they can form an important part of the analy- 
sis . The equation of state for air in equilibrium can be written 

P = p(h,p) (12) 

The transport properties are the total properties as defined by Hansen 
(ref. 9 ). Inspection of equations ( 9 ) and (10) reveals that the gas proper- 
ties occur in discrete, select groupings - that is, p, pp, and pp/cr. These 
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groupings have been correlated by Viegas and Hove (ref. 10) as functions of 
the static enthalpy, h, and the pressure, p, for equilibrium air. The corre- 
lations were fitted with analytical curves by Viegas and Hove for use on 
electronic digital computers. These curve fits vere used to obtain the solu- 
tions presented in the present report. Because the correlation curves of 
Viegas and Hove do not apply at temperatures belov 1000° K, t! perfect gas n 
relationships vere used for temperatures belov 1000° K. The "perfect gas” 
relations vere modified slightly to remove any discontinuity in properties 
at 1000° K. 


Solution of Equations 


Equations ( 9 ) and (10) are general for equilibrium flov except for the 
usual assumptions of a thin boundary layer and the disregard of radiative 
transport; but because the right sides of the equations are not similar, they 
are very difficult to solve and must be solved at many points along a given 
body shape. Hovever, for a cone or a flat plate vith isothermal vail, 
derivatives vith respect to £ are identically zero because neither the 
boundary conditions nor the flov properties vary in the £ direction (i.e., 
the solutions are similar) . Making use of this fact and setting 0 = 0 ve 
simplify equations ( 9 ) and (10) to 

(Zf Tin ) n + ff^ = 0 ( 13 ) 



Ug2 

He 



Zf T) f T|TJ 


= 0 


(14) 


Equations ( 13 ) and (l4) are coupled through the transport properties, 
state properties, and the function f and its derivatives. The coupling 
betveen the momentum equation and the energy equation for boundary-layer flov 
has been demonstrated to be veak or knovn under certain conditions. For 
example, reasonable stagnation-point heat-transfer rates (i.e., for |3 / 0) 
vere calculated by Lees (ref. 11) assuming the value of l to be unity. 

Also for a equal to unity the Crocco equation relates equations ( 13 ) and 
(l4) directly. This veak coupling and the Crocco equation have been used to 
obtain rapid solutions in the folloving manner. 


Assume a solution for equation (l4) . Equation ( 13 ) can then be con- 
sidered as a third-order ordinary differential equation vith knovn variable 
coefficients and its solution for given boundary conditions may be obtained 
readily (see next section for procedure used in solving the differential 
equations) . Similarly, equation (l4) can be solved for an assumed solution 
of equation ( 13 ) • Therefore, if an approximate solution to one of the 
equations is knovn, the exact solution may be obtained by iteration betveen 
the tvo equations. This process continues until the solutions converge; that 
is. 
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and 


( 15 ) 


i Si - Si_! i = <p a (16: 

where cp^ and cp 2 are required to approach zero with increasing number of 
iterations, i. Since it is impractical to check the functions f and g at 
many points, the derivatives of f and g were checked at the wall; that is, 

f mw i “ f mw i _ 1 = ^3 ^3 “*■ 0 (W' 


= ^ ^4 0 (18) 

i oo 

Because of the nonlinearities involved, this cannot he proved to he an 
absolutely convergent procedure, nor can it he shown to he unique. However, 
in cases where the solutions have heen compared with those obtained by other 
investigators using other procedures, the agreement is very good. 

There is still the problem of how to get the first approximate solution 
to start the iteration procedure. Consider equation (13); note that the only 
term which depends on equation (l4) is the ratio This is a function 

of the static enthalpy 

-Hi- = -BL (h) 

PwM-w Pw^w 

where h may he written as a function of g and the first derivative of f 


h = Big 


u e 2 - 2 

~ 1 


Now by the use of the Crocco equation (ref. 12) 


(l §w) 


h can be expressed as 


h = H e (l - g w ) f^ + H e g v + f ^ 2 

Equation ( 13 ) can now be solved independently of equation (l4) 


The process outlined above provides for very rapid convergence. A 
typical example of the convergence process is given in appendix A. 

Numerical analy sis proc edures employed.- Each of equations ( 13 ) and (l4) 
was solved by integration, using a fourth-order Adams-Moulton method (ref. 13 ) 
with variable step size and double precision arithmetic to minimize round-off 
errors. The procedure was to select an initial value of f^rj or grj at the 
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wall and integrate until t] became large , or some other condition such as 
f-q-q or g-n asymptotically approached zero or or g became excessively 

large . In any case a second value of f ^ or g T j was determined and the 
equation integrated again. This was repeated three times and the results 
linearly extrapolated to the outer boundary condition to obtain a fourth 
value. This procedure was repeated with smaller and smaller changes in the 
initial value until the solution had converged within prescribed limits. The 
limits prescribed for the results presented in this paper are 

. f ^ and g at large t] 1.0 ± 0.001 

f^ and g-^ at large r\ 0 ± 0.001 

f^ and g^ at T] = 0 ± 0.001 

These limits were found to give results which were accurate within ±1/2 per- 
cent. Large changes in these limits were found to have only a very small 
effect on the final results. See appendix A for more details on the effects 
of various quantities on the convergence. 


Heat-Transfer, Skin- Friction, and 
Boundary-Layer-Thickness Parameters 


The methods used to calculate heat transfer, skin friction, and various 
boundary-layer thicknesses from the solutions of equations ( 13 ) (l^-) are 
outlined in the following paragraphs. 


Heat transfer .- The heat transfer to a surface in equilibrium air in the 
absence of radiation is 


4 


V hi 
a-w “ 


St 

Sy 


w 


( 22 ) 


When the transformations in equation (7)j and the dimensionless variables 
(eq. (8)) are applied, equation (22) becomes 


k w Her n UePw 
" °p ^ 


( 23 ) 


This equation is general and can be used for either two- or three-dimensional 
flow and can be applied downstream of a slightly blunted nose or leading edge, 
if similarity is assumed to hold and 0 = 0. However, for the present purpose 
we consider only the case of a sharp cone with an isothermal wall. For this 
case p^, and u e are constant and the radius of cross section is given 

as r = x sin 0 C , where 9 C is the cone half-angle. With this restriction 
equation (23) may be written as 


9 



(214-) 


q. = 


k w He J3/2 


c„ x 
Pw 


5l lw 


where 


_ Pw%x 

■KW 

Mw 

This may be rearranged to form the Nusselt number, Nu, divided by the square 
root of the Reynolds number, 

Nu _ ^ c Py X _ ^w ^25 

^/R^ (He - H w )k w v 2 1 - g w 


All of the heat-transfer results presented herein are in this form. 

Skin friction .- The local shear stress, t, on the wall can he expressed 
as 


r 


8u 
^ 8y 


(26) 


When the transformations and dimensionless forms in equations (7) and (8) are 
applied, this becomes 


T 


^w^e 2 ** Pw 


fr nw 


(27) 


This equation is general in the same sense that equation (23) is for heat 
transfer. However, for the purpose of this paper we will again consider only 
the case of a sharp cone with isothermal wall. If properties consistent with 
such restrictions are used, and the shear stress is divided by the local 
dynamic pressure [(1/2) PgUg 2 ], the local skin- friction coefficient, c f , can be 
written as 


or 


where 


76 fT)1» 
f _ 7i;7T; 


Cfjsi - ^L£sii 


Re - Pe u e x /i Me 


(28) 


( 29 ) 
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and 


* 

= Pe^e/ Pw^w 


Note that different Reynolds numbers are used for heat transfer and skin 
friction to simplify the correlations. 


If we now integrate the local skin friction over the length of the cone 
and divide hy the surface area, we obtain an average skin- friction coefficient 

" 1 V6 f 

Cf feZ = ~ (30) 




•JTZ 


where 


Re 


x b 


P e u e x b 

^e 


Boundary-layer-thickness parameters There are several thickness 
parameters associated with boundary layers. The three most commonly used are 
the boundary-layer thickness, 6, defined as the distance, y, at which the 
boundary-layer-flow velocity has reached 0-995 the displacement thickness, 
5*, defined as 


5* 


-/ 


,y(u/ue = 0.995) 


1 - — 

P P 71 


ay 


(3i) 


and the momentum thickness, 0, defined as 


0 


^y(u/u e = 0 . 995 ) 





(32) 


These thicknesses for the case of cone flow, appropriately made dimensionless, 
were calculated by applying the transformations and dimensionless forms 
described earlier, as follows 


and 



(33) 


5* ^ 

X 



(34) 
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and 


#* 



where tj is the value of tj at which f = u/u e = 0-995* 
e *1 

The integrations indicated in equations (33) * (3^*)> and (35) were 
performed by use of the Simpson rule . 

Additional heat-transfer, skin- friction, and boundary-layer-thickness 
parameters are listed in appendix B, as well as formulas, for conversion from 
one form to another. All of the foregoing discussion on heat-transfer, skin- 
friction, and boundary-layer-thickness parameters is applicable to a flat 
plate merely by striking out the 3*s wherever they occur. 


RESULTS AND DISCUSSION 

Equations (13) and (l4) were solved for the following range of condi- 
tions: 

u^ 3*05 to 2h.k km/sec 

Ue 2 /2He 0 to 0-93 
p e 0.1 to 10 atm 

T w 1000° to 4200° K 

Solutions for the lower values of u^,, less than 10 km/ sec , were used to 
check the present solutions against others which included only dissociation 
effects. The upper limit was dictated by the maximum temperature for which 
the thermodynamic and transport properties of air are available. 

The parameter u e a /2E Q can be interpreted as a cone-angle parameter. 
Using impact theory for a cone (i.e., Ue = u^ cos 0 C ) and neglecting the small 
contribution of the free-stream static enthalpy to the total enthalpy, one 
can express this parameter as 

Ug 2 /2He = cos 2 9 C (36) 

where G c is the cone half-angle. (This must he modified for cones with 
small amounts of bluntness, but at the speeds of interest this modification 
is small.) By use of this relationship, values of u e z /2H e of 0, 0.25, 

0.50, 0.75, and 0.93 correspond to cone half-angles of 90 °, 60 °, 45°, 30°, 
and 15°, respectively. For the 90° cone, the shock wave will be detached and 
the assumption of conical flow will not apply. However, a different inter- 
pretation can be placed on the parameter in this case, since Ug 2 /2He = 0 
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also corresponds to solutions for zero Mach number with zero pressure gradient 
and very high free-stream static enthalpy. Such a condition may represent 
flow in the heat exchangers of a reactor , for example. 

Most of the solutions were obtained for a boundary-layer pressure of 
1 atmosphere. A few solutions were obtained at pressures of 0.1 and 10 
atmospheres to assess the effect of pressure. 

The wall temperatures were chosen to cover the range of probable interest 
for entry vehicles . 


Velocity and Enthalpy Profiles 


Typical examples of solutions to equations (13) and (ik-) are tabulated 
in table I. The variables of primary interest, u/ue, H/He, and h/h e are 
also plotted in figure 1 as functions of the dimensionless distance normal to 
the surface, t), for two cases. The velocity profiles are shown in figure l(a). 
These profiles are similar to those obtained by Van Driest (ref. l4) . The 
inflections exhibited by these curves result from transformation from y to T) 
and do not exist when plotted in the physical coordinates. A numerical 
inverse transformation was performed to verify this but it is too cumbersome 
to present here. 

The total-enthalpy profiles shown in figure l(b) are similar to those 
of reference l4. There is a slight overshoot in H/H e (i.e., g exceeds l) 

near the outer edge of the boundary layer for the larger values of Ue 2 /2He 

(difficult to detect at the scale plotted) . This overshoot is a result of a 
Prandtl number less than unity and the maximum temperature occurring near the 
center of the boundary layer. 

The static-enthalpy profiles are shown in figure l(c) . These profiles 
are indicative of the temperature profiles. Note that, as stated above, the 
maximum static enthalpy, or temperature, occurs near the center of the 
boundary layer for n e ^/2E e greater than about 0.5* These profiles are very 
similar to the temperature profiles of reference l4. Although it can not be 
detected in these figures, it is evident in table I that at the edge of the 

boundary layer, the static enthalpy has not returned entirely to the edge 

values for large Ue 2 /2He. This discrepancy is a result of the convergence 
criteria placed on u/ue and H/He^ which have to be combined to yield h/he* 
This does not, however, have any appreciable influence on the skin- friction, 
heat-transfer, or boundary-layer-thickness results. 


Comparison With the Crocco Equation 


The similarity of solutions of the momentum and energy equations for a 
Prandtl number of unity was first shown by L. Crocco for a perfect gas in 
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reference 12. The generalization of his result to a real gas in chemical and 
thermodynamic equilibrium follows in a similar fashion to yield equation (20) . 


A comparison of this result (eq. (20)) to the calculation for variable 
Prandtl number, the present solutions, is shown in figure 2. The close 
agreement probably accounts for the very rapid convergence of the calculation 
procedure employed, which uses the Crocco equation to help calculate the 
state properties and transport properties for the first approximate solution 
to the momentum equation. 

A commonly used way of extending the Crocco equation is to differentiate 
equation (20) with respect to the independent variable r\ to yield 

Srj = (l “ Sw) f (37) 

This equation provides a simple relationship between heat transfer and skin 
friction. However, it is not accurate for the present conditions, as can be 
seen from an examination of the initial slopes of the curves in figure 2. 

From analogy, this would imply that Reynolds analogy is questionable under 
these conditions. This was also noted in reference 8, 


Presentation and Correlation of Heat-Transfer Results 


The heat-transfer, skin- friction, and three boundary-layer-thickness 
parameters obtained from the solutions are listed in table II. Other informa- 
tion essential to devising correlations of these results is included. The 
effects of flow conditions and wall temperature on these results are discussed 
in the following paragraphs. Comparisons are made with other theoretical 
results where available. 

The heat-transfer results are shown in figure 3 in the form Nu/Jr^ as 
a function of velocity u^. The data in this form are easy to use, since the 
flow properties are evaluated at the wall where the effects of dissociation 
and ionization are least strong. The results are also correlated by this 
type of presentation. 

Correlation . - The general features of the correlation plots (fig. 3) are 
shown in sketch (b ) . Most of the present solutions were correlated by the 
straight line, Nu/Jr^ = K(p e ^T- w -)u O0 "‘°* 36 . Deviations from the correlation were 
noted with high wall temperature at low velocities and to a somewhat smaller 
extent with large cone angles at high velocities, as indicated by the dashed 
lines in sketch (b) . The reasons for these deviations will be discussed. 
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Sketch (h) 

Plotted in the manner of sketch (h) , the results are independent of the 
parameter 1 ^ 2 / 2 ^ (i.e., cone angle), except as noted. 

The effects of the wall temperature and pressure on the correlation 
enter through the parameter K(pe,T-w)j as illustrated in figure 4. Particu- 
larly at pressures above 1 atmosphere, the effect of pressure change is small 
compared to the effect of wall temperature change . This insensitivity to 
pressure above 1 atmosphere is a result of the small effect of pressure on 
the thermodynamic and transport properties, as may be observed in reference 10. 
In most cases of practical interest and at pressures of less than 1 atmos- 
phere the effects of nonequilibrium have to be considered, thereby introduc- 
ing additional complications not considered herein. It is thought that for 
most calculations the correlation parameter for a pressure of 1 atmosphere 
will give satisfactory results; this may be approximated by the equation 
K(p e ,T w ) = 0.280 T-^/*" 0,1335 , where the constant 0.280 has units of 
(km/ sec)° * 36 Ort 0 * 1335 . 

Deviations from the correlation .- There are two conditions under which 
significant deviations from the correlation curves occur. The first of these 
is when the wall enthalpy approaches the recovery enthalpy. This condition 
is characterized by low u^ and large u e 2 /2H e and is indicated by the 
dashed lines on the left side of sketch (b) . This deviation is a result of 
using the difference between stagnation enthalpy and wall enthalpy rather 
than the difference between recovery enthalpy and wall enthalpy in defining 
the Nusselt number. However, because of the simplicity of determining the 
stagnation enthalpy and the difficulties encountered in determining the 
recovering enthalpy, it was thought that the use of the stagnation enthalpy 
was more practical. 

The second condition under which deviations occur is when the air in the 
outer portion of the boundary layer becomes highly ionized. This happens only 
at speeds greater than 9 "to 10 km/sec and appears to be a strong function of 
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the parameter, ue 2 /2He; and pressure. The effect of the cone-angle parameter 
results from its effect on the static temperature and thus on the degree of 
ionization; decreasing ue 2 /2He (i.e., increasing cone angle), increases the 
static temperatures and thus increases the degree of ionization. This effect 
can he seen in figure 3(a) . The effect of pressure on this deviation from 
the correlation is one of changing the degree of ionization and may he seen 
hy comparing figures 3(b), 3(3-) > and 3(e). No appreciable effect of vail 
temperature on deviations due to ionization is evident, prohahly because the 
vail temperatures considered are relatively small in comparison to ionization 
temperatures . 

Comparison vith other theo retical results.- The present results vere 
compared vith those of references 7 and 8. For the purpose of comparison 
vith the results of reference 7 the logarithm of Nu/s/r^ vas plotted in 
figure 5(a) versus the logarithm of ( p e li e /pw^w) s t ( the subscr ipt st indi- 
cates that the ratio is evaluated at the stagnation condition behind a normal 
shock at the flight conditions in question.) The agreement of the present 
results vith the correlation curve of reference 7 is very good until ioniza- 
tion begins. At this point discrepancies betveen the present results and the 
extrapolated results of reference 7 appear and become very large vith increas- 
ing velocity. At a speed of 15 km/sec the present results are larger by a 
factor of 1.50 than the extrapolated results of reference 7. 

A comparison vith the theoretical results of reference 8 is shown in 
figure 5(b) . There the solid lines are the present correlation curves and the 
circles are results from reference 8. Note that for velocities of ^.36 and 
6.23 km/sec there is generally good agreement betveen the results. The agree- 
ment at 10.86 km/sec, however, is poor. This is probably a result of the 
transport properties used; in these particular results the calculated trans- 
port properties used (ref. 8) did not include the effects of ionization. 


Calculation of heat t ransfer to a cone.- Up to t his point the heat- 
transfer results have been presented in the form Nu//r^, vith all flow 
variables evaluated at either vail conditions or boundary-layer edge condi- 
tions. This form permits the results to be presented without restriction as 
to the solution of the inviscid cone-flov problem. In order to apply these 
results to a heat-transfer calculation, it is necessary to make some assump- 
tions about the inviscid cone flow. We will assume that all of the flow 
variables at the boundary-layer edge may be calculated by simple modified 
impact relations, where the pressure relationship is 


p e _ § c P m ax ftoU °° 2 sin2 + Poo 
Poo Poo 


(38) 


and the velocity relationship is 


u e = u oo cos G c 


(39) 
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where is the stagnation-point pressure coefficient; c p n]ax is "between 

1.8 and 1.9 for most cases. Using these relationships and the correlation 
equation shown in sketch (b) , we can write the local heat-transfer rate as 

* _ K(p e ,T w ) J PwMw 

If we now use an equation of state of the form 

p = p/zRT ( 4 l) 

and note that, for the wall temperatures and pressures considered, the 
compressibility factor z differs only slightly from 1 (we will therefore 
assume it to be identically l), that H e » and that p e » p^, we can 
rewrite equation ( 4 o), using the expression K(T w ,pe) = 0.280 Tw 0 ' 1335 ^ as 



where the constant 0.268 has the units (km/sec) ‘ 36 °K ~° ,X335 # The terms in 
the first parentheses represent the dependence on wall temperature (the 
effect of pressure being very small on any of these terms) . This group of 
terms is nearly constant, varying only a few percent for wall temperatures 
from 1000 ° to 4200 ° K. The term in the second set of parentheses is the 
altitude dependence; the third group, the geometry dependence; and the last 
group, the velocity dependence. 

In practice, the tips of the cones will generally be blunt either by 
design or as a result of ablation. The influence of small amounts of tip 
blunting on the skin friction and heat transfer is considered in appendix C. 


y C ° S X — Uoo° ■ 14 (H e - H w ) (b0 ) 


Presentation and Correlation of Skin-Friction Results 


Typical skin-friction results are shown in figure 6. Formulas for 
conversions to other parametric forms are given in appendix B. 

Correlation .- The skin- friction results were not so readily correlated 
as the heat-transfer results. A useful correlation was obtained when the 
parameter Cf^yRe was plotted versus u e 2 / 2 H e as shown in figure 6. Note 
that Cf^/R^ decreases with increasing u e 2 /2H e (i.e., with decreasing cone 
angle) . The equation of the correlation curves shown, which fits the results 
fairly well in most cases, is 

c f -/Re = a - b(ue 2 /2H e ) m (^ 3 ) 
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where a is a function of the wall temperature as shown in figure 1, b is 
a constant equal to 0.506, and m is a constant equal to 2.2. This simple 
correlation curve represents the calculated values within ±10 percent in most 
cases . 

Deviations from the correlation .- The main deviations from the correlation 
curve (fig. 6) occur at the same conditions at which deviations from the 
heat-transfer correlation occur, namely, (l) near-zero heat transfer and (2) 
when a large fraction of the outer part of the "boundary layer becomes appreci- 
ably ionized. The first of these effects is evident in figure 6(b) . Here the 
wall enthalpy has become a large fraction of the stagnation enthalpy (i.e., 
nearing zero heat transfer) for the 4.56 km/sec results, causing them to 
deviate from the 9-12 km/sec results; this is also evident in results for 
higher wall temperatures . 

The effect of ionization is evident to some degree at all wall tempera- 
tures and pressures considered. It results in an increase in the skin- 
friction parameter with decreasing ue 2 /2He for the higher speeds - that is, 
higher total enthalpies . Also the lower pressure results are not so well 
correlated as the higher pressure results, an effect which is again attribut- 
able to ionization. 

Comparison with other theoretical results.- Comparison with the theo- 
retical perfect-gas results of VanTllriest (ref. l4) is shown in figure 8(a). 
The curve from reference l4 is for a ratio of T w /Te (or h w /h e ) of 0.25* 

The present results agree fairly well with those of reference l4. The slight 
increase in the skin-friction parameter for the present results is attribut- 
able to real-gas effects. The discrepancy is greater at the smaller values 
of Ue 2 /2He because of higher temperatures in the boundary layer. The 
inflection in the skin-friction curve near 1 u e 2 /2E e equal to zero is probably 
an effect of the beginning of ionization. 

Comparison with the results of reference 8 and the "reference enthalpy" 
method of reference 15 is shown in figure 8(b) . The present results agree 
well with both of these sets of results for nearly the entire range of con- 
ditions considered in the present report. The exception was for conditions 
when the boundary-layer air was appreciably ionized. 

It would appear that the reference enthalpy method would predict skin 
friction more accurately than the present correlation equation, but would 
probably be more difficult to apply. 

Calculation of the friction drag coefficient on a cone.- The skin- 
friction drag coefficient, based on the base area, is obtained when the local 
skin friction is integrated over the cone as follows: 


c D f 


2 sin 6 C cos G c 
rb 2 Poo"Uoo 2 



m 
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r 


Substituting equation (39) 


into equation (44) yields 




2 sin 9c cos 3 Q jg p e 
r^ 2 Poo 



(45) 


No w using the correlation equation for the skin-friction parameter and equa- 
tion ( 36 ), integrating and rearranging , we obtain 


where 


c D f 


4 cos 5/2 0 C ( a 
3 sin 0 C 


0.506 cos 4,4 0 C ) / PeM-e 

V^ooXb V Poo,io0 


Roo^b 


PpO^h 

Moo 


(46) 


It is very difficult to ascertain the separate effects of various conditions 
on Qj) because of the dependence of p e p e on both the free-stream velocity 
and the cone angle . 

Some general effects of small amounts of bluntness are discussed briefly 
in appendix G . 


Boundary-Layer-Thickness Parameters 


Three boundary-layer-thickness parameters, the boundary-layer thickness 
6 , the displacement thickness 5*, and the momentum thickness 0 , made dimen- 
sionless by (l/xX/5^ were calculated for various values of Ue 2 / 2 He* 

Boundary-layer thickne s' s . - The boundary-layer thickness out to 
u/ue = 0.995 is shown in figure 9 . For a given local Reynolds number the 
boundary layer increases in thickness with increasing u^ 2 / 2 He (i.e., decreas- 
ing cone angle). This general trend holds for all velocities, wall tempera- 
tures, and pressure levels considered. The effect of pressure level was 
small. The effects of wall temperature and free-stream velocity on the 
boundary-layer thickness appear to be closely related and their combined 
effect appears to correlate with g w (i.e., 2h w /u O0 2 ) . This is shown in 
figure 10. Here some typical values for a dimensionless boundary-layer thick- 
ness are plotted versus g w . Except for some scatter at small values of g v , 
increasing g w increases the boundary-layer thickness. 

The effect of ionization appears to be one of thickening the boundary 
layer; this occurs at high speeds and for small values of u e 2 /2H e . Note in 
figure 9 the boundary-layer thickness for 15*20 km/sec increases with decreas- 
ing UgS/^He below u^/^He = 0-5- 
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A "brief comparison of the "boundary-layer thickness with the theoretical 
results of reference l4 is shown in figure 9* The curve from reference l4 is 
for a ratio of t%/he of 0.25* No attempt at putting the present results 
into a constant hv/h e form was made because of the small variations in the 
boundary-layer thickness with the various parameters . The agreement between 
the results of reference l4 and the present results is good. The agreement 
would be better if the effect of changing h w /h e were considered. 

Displacement thickne ss.- Typical results of boundary-layer displacement 
thickness are plotted versus Ue 2 /2He in figure 11. The general trends of 
the results, with the various parameters u^/^He, T v , Uoo, and p e are very 
similar to those exhibited by the boundary-layer thickness, except the effect 
of ionization is not noticeable on the displacement thickness. There is one 
very interesting fact evident in this figure - that is, the displacement 
thickness appears to get very small and possibly take on even negative values 
for small values of the parameter Ue 2 /2H e when the value of g w is very 
small. This phenomenon may be- of interest in the design of high-speed inlet 
systems or forced gas heat exchangers in nuclear reactors . Under these con- 
ditions the mass flow through the system might be larger than one would 
expect . 

Momentum thickness .- Typical results for boundary-layer momentum thick- 
ness are plotted in the same manner as the other boundary-layer thicknesses 
in figure 12. The general trends exhibited by these results are very similar 
to the trends exhibited by the skin- friction results. The reason for this 
similarity is that the momentum thickness is, by definition, proportional to 
a momentum defect (or loss) from the main flow over the cone. This loss has 
to appear as drag to the surface over which the boundary layer flows . This 
fact has been used for years in subsonic wind tunnels to determine drag. 


SUMMARY OF RESULTS 


The boundary-layer equations in axisymmetric flow with zero pressure 
gradient for air in chemical and thermodynamic equilibrium have been solved 
for a wide range of conditions. From these solutions the heat transfer, skin 
friction, boundary-layer thickness, displacement thickness, and momentum 
thickness for cones in hypervelocity flight have been computed. The follow- 
ing are some of the results of these calculations. 

Good correlation of the heat-transfer results was obtained when the 
Nusselt number divided by the square root of the Reynolds number, R w , was 
plotted versus the free-stream velocity. The effect of the cone-angle param- 
eter on this correlation was negligible. The effects of wall temperature and 
pressure level were accounted for in the correlation. Two regions of poor 
correlation were noted, one near zero heat transfer and the other at condi- 
tions where the degree of ionization was very high. The present results 
indicate heat-transfer rates which are higher than extrapolations of earlier 
theoretical results which include only dissociation and neglect ionization. 
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A fair correlation was found for the skin- friction results. This corre- 
lation allows for the simple estimation of the friction drag. A letter 
method for calculating skin friction is the reference enthalpy method; how- 
ever, it is somewhat more difficult to use. There were no major effects of 
dissociation and only small amounts of ionization on the skin-friction param- 
eter, c f^/Re". However, large amounts of ionization tend to increase the value 

of Cfv fe. 

No attempt was made at correlating the houndary-layer-thickness param- 
eters; however, from comparisons with solutions for a perfect gas, it was 
noted that there were only minor effects of dissociation and ionization. One 
interesting point evident in the displacement-thickness results was that very 
small and even negative values were indicated under conditions which simulated 
subsonic flow at very high static enthalpies with cold walls. These results 
may he of some interest in the design of hypersonic inlets and nuclear-reactor 
heat exchangers . 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., June 15 , 196^ 


21 



APPENDIX A 


CONVERGENCE OF THE ITERATION PROCEDURE 


The over-all convergence procedure of iterating "back and forth between 
the momentum equation and the energy equation is illustrated in figure 13* 
Values of fq and g at large t) are plotted against corresponding values of 
f TVOw an( ^ Stjw’ *^ ie ^- aslie( ^ li nes indicate the selected convergence limits. The 
iteration procedure consists in obtaining a solution to the momentum equation 
using the Crocco equation to determine the state properties and transport 
properties. This is indicated by I in figure 13(a). We then proceed to 
obtain solutions for three arbitrary values of g-q and plot the value of g 
obtained for large r\ versus the guessed value of gq - This is indicated by 

II in figure 13(b) . We then draw the line connecting these results and 
extrapolate toward the correct outer boundary condition to obtain a fourth 
value for gq v and solve the energy equation (point designated k g) . We then 
return to the momentum equation and perform a similar operation but now 
replace the Crocco assumption for calculating the transport properties and 
state properties with the results from the previous solution to the energy 
equation. This process continues, iterating back and forth between the momen- 
tum-equation and the energy equation, until both solutions have converged 
within the prescribed limits . It should be noted that the high degree of 
accuracy obtained by the linear extrapolation was not anticipated, and is 
probably a result of the good approximation of the initial solution to the 
momentum equation from the Crocco equation. 

EFFECTS OF STARTING VALUES ON THE CONVERGENCE PROCEDURE 

There was no effect, at least to four significant figures, of starting 
values on the final converged solutions. This is illustrated below. 

EFFECT OF STARTING VALUES ON CONVERGENCE 

Voo = 30,000 ft/sec T w = 500° K 
p e = 1 atm Ue = 0 


Starting values 

Converged values 

frjT](0) 

Stj(O) 

f qq( 0 ) 

Stj(0) 


s n (°°) 

0.320 

0.210 

( 9 ) 

0.3225 

(11) 

0.2135 

0 . 4487-4 

O.IIO8-3 

• 300 

.215 

( 9 ) 

.3226 

(11) 

.2135 

. 4487-4 

.1108-3 

O 

00 

OJ 

.220 

( 9 ) 

.3226 

(11) 

.2135 

. 4486-4 

r 

.1x07-3 

0 

D'— 

OJ 

.225 

(10) 

.3226 

(11) 

.2135 

. 4486-4 

.1108-3 

.260 

.230 

(10) 

.3226 

(12) 

.2135 

. 4487-4 

.1108-3 
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The numbers in parentheses indicate the total number of iterations for the 
given equation. Note that there is not much effect of starting values on the 
number of iterations; this is again a result of the high degree of linearity 
exhibited by a plot of f^ or g^ at large values of t) versus f^ or g^ 
at the wall . 


EFFECT OF THE SIZE OF THE CONVERGENCE CRITERIA ON THE SOLUTIONS 


The effect of the size of the convergence criteria on the converged 
solutions was very small. This is illustrated by a typical example below. 


Convergence 

criteria 


grj(0) 


St](°°) 

lxicr 3 

0.424635 

0.286245 

0.1721X10- 3 

0.2299X10" 3 

5X10 -4 

. 42 I 1-592 

.286245 

•5022X10' 5 

.1246XL0" 3 

2X10“ 4 

.424592 

. 285902 

.5807X10" 6 

•4377X10" 4 


Note- that as the convergence criteria on u/ue and H/H e at the outer edge of 
the boundary layer is decreased from 1 ±1X10“ 3 to 1 ±2X10~ 4 , the converged 
values of frjr^ change only in the fifth significant figure and g^ 
changes only in the fourth significant figure. There is some improvement in 
these quantities at large values of r|; however, they are all within the pre- 
scribed limits of 0 ±1X10~ 3 . This insensitivity to the size of the conver- 
gence criteria can again be attributed to the high degree of linearity 
exhibited by the plots in figure 13 . 

The process outlined above provides for very rapid convergence; solu- 
tions were obtained in less than one-half minute on an IBM 7090 electronic 
digital computer. The process has also been used successfully for stagnation- 
point solutions, for cases involving both moderate pressure gradients and a 
finite velocity at the edge of the boundary layer and cases involving moder- 
ate amounts of blowing or suction at the wall. None of these results are con- 
tained in this report. 
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APPENDIX B 


HEAT-TRANSFER AND SKIN-ERICTION PARAMETERS 
OF COMMON INTEREST 


There are other heat-transfer and skin-friction parameters in vise Besides 
those described in the text. Some of these will be described below, special- 
ized to the case of cone flow. 


The Nusselt number divided by the square root of the Reynolds number, 
R w , as described in the text, is 



(Bl) 


This may he rewritten in terms of the local Reynolds number, R e , as 


where 


Nu = / 3 ( PePe\ Pw g l v 

v 2 WW Pe 1 - gw 


Re = 


Pe 


(B2) 


Conversion from the parameter described by equation (Bl) to that described by 
equation (B2) is accomplished as follows 


Nu 


Nu 


L W 


7 V2 

e P^ 


(B3) 


The Stanton number, St, is another commonly used heat-transfer parameter 
and is defined as 


St = 


(H e _ H w ) pe^e 


(B4) 


The Stanton number usually appears multiplied by the square root of a Reynolds 
number as 


St 



1 - g-w 


- 1/2 

^e 


-l 

<*w 


(B5) 


2b 



which can he put in terms of the parameter used in the present report as 


‘ (jfe) V 1 <B6) 

The quantity appearing in parentheses in the foregoing equation is the 
heat-transfer parameter used in the present report; all other quantities 
appearing in the above equations are tabulated in table II. 

The only variation in the skin- friction parameter, c^/r^, which may be 
of interest, is to substitute R w for Re as follows 

Cf^fW = (c f VR^) I e l/s ^ (B7) 

where the term in parentheses is the parameter as used in the present report. 
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APPENDIX C 

EFFECTS OF SMALL AMOUNTS OF NOSE BLUNTING 
ON SKIN FRICTION AND HEAT TRANSFER 


In the strictest sense nose blunting completely invalidates the similar- 
ity assumption used in obtaining the solutions presented in this report; 
however, for small amounts of blunting and sufficiently far downstream from 
the nose (i.e., downstream of the large pressure gradient), the assumption of 
similarity again will probably be good. We still have to modify the pure 
cone results in two ways. First, the value of the parameter u^/Hq has 
been altered because the flow has passed through a normal shock resulting in 
somewhat lower edge velocities. Second, the transformed coordinate £ can 
no longer be obtained directly because the integral used in defining £ 
cannot be integrated without specification of the nose geometry which will 
dictate the pressure history along the body and the flight conditions. 
Therefore, the heat transfer and skin friction would have to be obtained by 
applying equations (23) and (27) of the text, being sure that the values of 
g^ and fr)^ are selected for the proper values of u e 2 /H e . 

The above discussion has disregarded the problem of the vorticity inter- 
action between the shock layer and the boundary layer. Because of the involved 
nature of the vorticity interaction, it is not considered here. It could, 
however, influence the results to a degree, depending on the amount of blunt- 
ing, the shock-layer thickness, and the boundary-layer thickness. 
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TABLE I.- TYPICAL SOLUTIONS 


(a) Uco = 15*2 km/sec; p e = 1 atm; T w = 1000 ° K; Ue 2 / 2 He = 0.25 


n 

’ f 


f nn 

n 

S 

Sti 

h/h e 

0 

0 

0 

0.246 

0 

0.009 

0.226 

0.012 

.008 

.000 

.002 

.252 

.016 

.013 

.243 

.017 

.055 

.000 

.015 

•299 

.062 

.025 

.289 

• 033 

.094 

.001 

.027 

• 352 

.i 4 i 

.051 

•359 

.068 

•133 

.003 

.042 

. 4 i 4 

.266 

.102 

.464 

.135 

.172 

.005 

.060 

.473 

.344 

.i 4 i 

• 534 

.187 

.211 

.007 

• 079 

.523 

.422 

.186 

.626 

.246 

.250 

.011 

.100 

.562 

.547 

.278 

.868 

.367 

.289 

.015 

.123 

• 593 

.609 

.337 

1.024 

.444 

.328 

.020 

.147 

.620 

.664 

.394 

• 993 

.518 

.391 

.031 

.187 

.658 

.703 

.429 

.801 

.560 

.469 

.047 

.239 

.696 

• 730 

.449 

.655 

.580 

.547 

.068 

.295 

.738 

.746 

.459 

.620 

.583 

.617 

.091 

.349 

.781 

.769 

.473 

• 593 

• 590 

.656 

.105 

.380 

.806 

.801 

.491 

• 579 

.607 

.695 

.121 

.412 1 

.825 

.832 

.509 

• 575 

.622 

• 730 

.136 

. 44 i 

.833 

.895 

.545 

• 575 

.678 

• 750 

.144 

• 457 

.834 

• 957 

.581 

• 576 

• 705 

.785 

.161 

.487 

.833 

r.066 

.644 

.568 

.780 

• 793 

.165 

.493 

.832 

1.191 

• 713 

.538 

1 .870 

.824 

.181 

• 519 

.827 

1.316 

• 777 

.482 

•946 

.895 

.219 

• 577 

.809 

1 . 44-1 

• 833 

, 4 o 8 

1.000 

• 973 

.267 

• 639 

.776 

1.566 

.879 

.326 

1.036 

1.051 

.319 

.698 

• 73 Q 

1.691 

• 915 

.248 

1.058 

1.128 

.376 

• 752 

.671 

1.848 

.947 

.167 

1.074 

1.286 

.501 

.846 

• 523 

2.004 

.968 

.108 

1.052 

1.629 

.815 

.965 

.185 

2.160 

.981 

.067 

.998 

1.941 

1.123 

•995 

.034 

2 . 348 

• 991 

.036 

.991 

2.254 

1.435 

•999 

.003 

2.723 

.998 

.009 

.998 

2 .566 

1.747 

1.000 

.000 

2.973 

1.000 

.003 

1.000 

3-223 

2.403 

1.000 

.000 

3.254 

3-535 

1.000 

1.000 

.001 

.000 

1.000 

1.000 
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TABLE I.- TYPICAL SOLUTIONS - Continued 


(t>) Uoo = 15*2 km/sec; p g = 1 atm; T w = 1000° K; u e 2 /2H e = 0.5 


4 

- f * 

f 4 

f 

44 

4 

g 


h/h e 

0 

0 

0 

0.260 

0 

0.009 

0.213 

0.018 

.016 

.000 

.004 

• 273 

.016 

.013 

.227 

.026 

.086 

.001 

.026 

•351 

.062 

.024 

.269 

.047 

.133 

.003 

.044 

.422 

.125 

.043 

.321 

.084 

.180 

.005 

.065 

.492 

.219 

.076 

■ 391 

.148 

.219 

.008 

.086 

.541 

.312 

.117 

.471 

.227 

.258 

.012 

.108 

.580 

. 4 o 6 

.165 

•550 

.318 

•297 

.017 

.131 

.612 

.500 

.221 

.643 

.425 

• 336 

.022 

.155 

.638 

.594 

.286 

• 753 

• 548 ! 

• 375 

.029 

.181 

• 659 

.68 7 

.362 

.873 

.691 

• 453 

.045 

.234 

.696 

•734 

. 4 o 4 

.918 

-753 

• 531 

.065 

.289 

.726 

.781 

.448 

.928 

.812 

.609 

.090 

.347 

• 751 

.828 

.491 

.8 86 

.862 

.687 

.120 

. 4 o 6 

• 773 

.875 

• 530 

.811 

.895 

.766 

.154 

.468 

•790 

.922 

.567 

• 735 

.915 

.844 

•193 

.529 

•793 

.961 

.594 

.685 

.908 

.922 

.236 

• 591 

• 774 

.992 

.615 

.666 

.881 

1.023 

• 300 

.66 7 

.719 

1.039 

.646 

.645 

.847 

1.102 

•355 

.721 

.663 

1.102 

.686 

.618 

.852 

1.211 

.437 

.788 

• 571 

1.164 

• 723 

.588 

.825 

1.367 

.567 

.866 

.427 

1.289 

.792 

.511 

.834 

1.523 

•707 

.922 

.291 

l. 4 i 4 

.850 

.418 

.850 

1.680 

.854 

• 959 

•179 

1.539 

.897 

• 321 

.874 

1.836 

1.006 

.980 

• 099 

1.664 

• 931 

• 233 

.902 

1.992 

1.60 

•991 

.049 

1.789 

• 955 

.160 

.928 

2.148 

1.315 

•996 

.022 

1.914 

• 972 

.105 

• 952 

2.305 

1.471 

• 999 

.009 

2.102 

.986 

• 053 

• 974 

2.617 

1.783 

1.000 

.001 

2 . 4 l 4 

.996 

.017 

.992 

2.930 

2.096 

1.000 

.000 

2.727 

• 999 

.006 

.998 

3-242 

2.408 

1.000 

.000 

3-039 

1.000 

.002 

1.000 

3.617 

2.783 

l'.OOO 

.000 

3.602 

1.000 

.000 

1.000 
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TABLE I.- TYPICAL SOLUTIONS - Continued 


(c) Uoo = 15*2 km/sec; p e = 1 atm; T w = 1000° K; Ug 2 /2H e = 0-75 


n 

f 

f 4 

f T)Tl 

n 

g 

Sr, 

h/h e 

0 

0 

0 

0.274 

0 

0.009 

0.214 

0.037 

.016 

.000 

.064 

.287 

.016 

.013 

.229 

.052 

.070 

.001 

.022 

.349 

.047 

.020 

.258 

•079 

.109 

.002 

.036 

. 4 o 6 

.094 

•033 

•299 

.128 

.148 

.oo 4 

.053 

.468 

.172 

•059 

.367 

.228 

.187 

.006 

.073 

.526 

.266 

.098 

.449 

• 392 

. 22.1 

.009 

.094 

.574 

• 359 

.143 

.528 

.545 

.266 

.013 

.118 

.613 

• 453 

• 197 

.608 

.746 

• 312 

.020 

.147 

.648 

.547 

.257 

.685 

.963 

•391 

.033 

.200 

.692 

. 64 i 

• 325 

.752 

1.180 

.469 

.051 

.255 

• 723 

.734 

.398 

.798 

1-397 

.547 

.731 

.312 

.745 

.828 

• 473 

.813 

1.600 

.625 

•998 

.371 

• 759 

.922 

.549 

.19k 

1.783 

• 703 

.131 

.431 

• 763 

1.016 

.621 

.746 

1.927 

.781 

.167 

.490 

• 757 

1.109 

.688 

.678 

2.032 

• 859 

.208 

.549 

• 74 o 

1.172 

.729 

.626 

2.012 

• 937 

.253 

.606 

• 711 

1.234 

.766 

• 572 

I.962 

1.016 

.302 

.660 

.670 

1.297 

.800 

.516 

I.893 

1.094 

.356 

.710 

.620 

1-359 

.831 

.461 

1.812 

1.172 

.413 

.756 

.562 

1.422 

.858 

. 4 o 8 

1.717 

1.250 

.474 

.798 

.501 

1.484 

.882 

• 357 

1.618 

1.328 

.538 

.834 

.438 

1.609 

.921 

.265 

1-597 

l. 4 o 6 

. 6 o 4 

.866 

. 377 

1.766 

• 954 

.172 

1.566 

1.547 

■ 729 

.912 

.276 

1-953 

•979 

• 093 

1.421 

1-703 

.875 

.947 

.183 

2 .l 4 l 

•991 

.045 

1.274 

1.859 

1.025 

• 970 

.114 

2.328 

•997 

.019 

1.165 

2.016 

1.177 

.984 

.066 

2.516 

•999 

.007 

1.091 

2.172 

1-332 

• 992 

.036 

2.703 

1.000 

.002 

l.o 48 

2.328 

1.487 

•996 

.018 

2.891 

1.001 

.001 

1.028 

2.516 

1.674 

•998 

.007 

3.078 

1.001 

.000 

1.016 

2.828 

3-203 

1.987 

2.361 

•999 

1.000 

.001 

.000 

3.234 

1.001 

.000 

1.010 



TABLE I.- TYPICAL SOLUTION'S - Continued 
(d) Ucq = 15*2 k m/sec; p e = 1 atm; T w = 1000° K; u e 2 /2H e = 0-93 


n 

f 

f n 

f nn 

0 

0 

0 

0.286 

.008 

.000 

.002 

.293 

• 055 

.000 

.017 

.344 

.102 

.002 

.035 

.412 

.148 

.oo4 

.056 

.490 

•195 

.007 

.081 

.560 

.242 

.011 

.108 

.6l4 

.289 

.017 

.138 

• 655 

• 375 

.031 

.197 

• 707 

.469 

.053 

.265 

.742 

.562 

.081 

• 335 

.760 

.656 

.116 

.407 

.763 

• 750 

.157 

.478 

• 750 

.844 

.205 

.547 

.720 

• 937 

.260 

.612 

.676 

1.031 

.320 

.673 

.620 

1.125 

.386 

.728 

• 556 

1.219 

.457 

■ 111 

.489 

1.312 

.532 

.820 

.421 

1.4o6 

.610 

.856 

.356 

1.500 

.692 

.887 

.296 

1-594 

.776 

.912 

.242 

1.687 

.863 

• 932 

.195 

1.812 

.981 

• 953 

.142 

2.000 

1.162 

.974 

.084 

2.187 

1.346 

.986 

.047 

2-375 

1.531 

.993 

.025 

2.562 

1.718 

■ 991 

.013 

2.750 

1.905 

• 99 8 

.006 

3-125 

2.279 

■999 

.001 

3-500 

2.654 

1.000 

.000 

4.062 

3.217 

1.000 

.000 


7 

g 

6t) 

h/h e 

0 

0.009 

0.217 

0.131 

.016 

.013 

.232 

.186 

.062 

.025 

•277 

• 353 

.109 

• 039 

• 320 

• 541 

.l4i 

.049 

•379 

.658 

•234 

.086 

• 439 

l.l4i 

.328 

• 131 

• 525 

1.716 

.422 

.184 

.606 

2.375 

.516 

.245 

•679 

2.983 

.609 

• 311 

• 735 

3-508 

• 703 

.382 

.768 

3.965 

• 797 

.454 

• 774 

4.283 

.891 

.526 

• 756 

4.477 

.984 

• 595 

• 719 

4.523 

1.078 

.661 

.668 

4.465 

1.172 

.720 

.607 

4.267 

1.266 

• 774 

• 539 

4.015 

1.391 

.836 

.444 

3.921 

1.578 

.906 

.306 

4.008 

1.766 

• 952 

.190 

3.864 

1.828 

.963 

.158 

3-303 

2.016 

.985 

.084 

3-020 

2.203 

.996 

•039 

2.688 

2.391 

1.001 

.015 

2.233 

2.578 

1.002 

.003 

1.710 

2.766 

1.002 

-.002 

1.398 

2.953 

1.002 

-.003 

1.214 

3*l4l 

1.001 

-.002 

1.094 

3.328 

1.001 

-.002 

1.067 

3-578 

1.001 

-.001 

1.04-1 

4.016 

1.000 

-.000 

1.000 
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TABLE I.- TYPICAL SOLUTIONS - Continued 


(e) Ucq = 15-2 km/sec; p e = 1 atm; T w = 3000° K; Ue 2 / 2H e = 0.5 


\ 

f 

... f 4 


0 

0 

0 

0.308 

.008 

.000 

.002 

.315 

.031 

.000 

.010 

• 34o 

.070 

.001 

.024 

.380 

.117 

.002 

.043 

.423 

•195 

.007 

.078 

• 479 

.273 

.015 

.117 

.518 

• 367 

.028 

.168 

•553 

.445 

.043 

..212 

.578 

• 555 

.070 

•277 

.607 

•Til 

.120 

• 375 

.643 

.852 

.180 

.467 

.666 

• 930 

.218 

• 519 

.666 

1.008 

.261 

• 570 

• 653 

1.055 

.288 

.601 

.639 

1.109 

.322 

.635 

.619 

1.203 

.384 

.691 

• 577 

1-359 

•499 

• 775 

.490 

1.516 

.625 

.844 

-392 

I.672 

.762 

.897 

.294 

1.828 

.905 

• 936 

.206 

1.984 

1.054 

.963 

• 135 

2 .l4l 

1.205 

•979 

.082 

2.297 

1.359 

.989 

.o46 

2.453 

1.514 

.994 

.024 

2.641 

1.701 

• 997 

• 010 

2-953 

2.013 

•999 

.002 

3.266 

2.325 

•999 

.000 

3.578 

2.638 

•999 

.000 

3.891 

2.950 

•999 

.000 

4.266 

3.324 

•999 

.000 


4 

g 

gT] 

h/h e 

0 

0.035 

0.239 

0.070 

.016 

.039 

.247 

.078 

.047 

.047 

.263 

.094 

.094 

.060 

.287 

•119 

.187 

.089 

• 335 

.176 

.281 

.123 

• 384 

.234 

.4o6 

•175 

.457 

• 336 

•531 

.238 

.547 

.448 

.656 

• 313 

.656 

.581 

•766 I 

• 390 

.744 

• 703 

.828 

• 437 

■ 753 

• 733 

.891 

.483 

• 711 

.748 

• 953 

•525 

.641 

.781 

1.016 

.563 

• 577 

.801 

1.039 

• 577 

.561 

• 793 

1.070 

•594 

.549 

.785 

1.125 

.624 

• 535 

• 771 

1.187 

.657 

-521 

• 713 

1-297 

.712 

.492 

.712 

1.422 

• 771 

.445 

• 737 

1.516 

.810 

.4oi 

.744 

1-703 

■ 877 

• 303 

.827 

2.016 

.947 

.154 

• 936 

2.328 

•979 

.064 

.980 

2.578 

.990 

.029 

■992 

2.828 

•995 

.013 

.996 

3-078 

.998 

.006 

.998 

3.328 

•999 

.003 

• 997 

3.578 

•999 

.001 

•997 

3.891 

1.000 

.001 

1.002 

4.266 

1.000 

.000 

1.002 


33 


TABLE I.- TYPICAL SOLUTIONS - Concluded 


(f) Uoo = 24.32 tan/sec; p e = 1 atm; T w = 1000° K; ue a /2H e = O .93 


n 

f 

f n 

frjT] 

n 

g 

g n 

hAe 

0 

0 

0 

0.239 

0 

o.oo 4 

0.181 

0.051 

.004 

.000 

.001 

.245 

.016 

.007 

.208 

-099 

.023 

.000 

.006 

.283 

.078 

.023 

.305 

.328 

.043 

.000 

.012 

.334 

•l 4 i 

.045 

.396 

. 64 i 

.066 

.001 

.021 

. 4 o 6 

.203 

.073 

.502 

1.037 

.102 

.002 

.037 

.502 

.266 

.108 

.659 

1.524 

.129 

.003 

.051 

.554 

• 312 

.143 

.818 

2.007 

.168 

.005 

.074 

.608 

• 336 

.163 

.869 

2.255 

.230 

.011 

.114 

.676 

. 4 o 6 

.218 

.667 

2.940 

.309 

.022 

.170 

• 767 

.465 

•257 

■ 699 

3.286 

• 371 

.034 

.221 

.852 

• 590 

.353 

.833 

4.392 

.400 

.o 4 i 

.246 

• 8 77 

.887 

.623 

.906 

8.093 

.420 

.o 46 

.264 

.889 

1.184 

.850 

.584 

11.212 

.439 

.052 

.281 

.898 

1.309 

.912 

• 375 

11.975 

.459 

.057 

•299 

.906 

1.363 

.929 

.259 

12.079 

.494 

.068 

• 331 

.918 

l. 44 i 

.945 

.174 

12 . o 4 o 

.533 

.082 

.367 

.926 

1.520 

• 957 

.143 

11.877 

.572 

•097 

.403 

.929 

1.598 

.968 

.120 

ll . 666 

.611 

.113 

.439 

.927 

1.676 

• 976 

.097 

11.378 

.68 2 

.147 

.504 

■909 

1.754 

.983 

.076 

10.664 

.760 

.189 

.574 

.866 

1.832 

.988 

.058 

9-733 

.838 

.236 

.639 

• 8 o 4 

1.918 

.992 

.042 

8.743 

.916 

.288 

.699 

.728 

1.996 

■ 999 

.030 

7.720 

.994 

.345 

• 752 

.643 

2.074 

■ 997 

.021 

6.727 

1.072 

. 4 o 6 

• 799 

.556 

2.223 

■999 

.009 

5.788 

1.150 

.470 

.839 

.471 

2.379 

1.000 

.003 

4.932 

1.260 

•564 

.885 

.362 

2.535 

1.000 

.000 

3.879 

1.353 

.649 

.915 

•279 

2.691 

1.000 

-.001 

3.162 

1.510 

• 795 

• 950 

.172 

2.801 

1.000 

-.001 

2.295 

1.744 

1.021 

.978 

.082 

3.098 

1.000 

-.001 

1.578 

2.057 

2.807 

1.330 

2.078 

.994 

1.000 

.027 

.001 

3-223 

1.000 

.000 

1.159 
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table II.- TABULATION OF FLOW CONDITIONS AND PRIMARY RESULTS 
(a) p e = 1 atm; ^=1.06 km 2 /sec 2 ; T w = 1000° K; a w = 0-71 


"Uco; 

km/ sec 

He, 

km 2 / sec 2 

•Ue 2 /2He 

g w Xl0 2 

he, 

km 2 / sec 2 

h w 

— X 10 
he 


P e 

— X 10 
Pw 

h e 

S-n 

V 

f Wv 

Nu/n/r^ 

cfs/Re 

1 ^ 


1-s 

4.56 

10.39 

0.0 

10.20 

10. 40 

1.020 

0.44l 

1.602 

0.6l 

0.301 

0.371 

o. 4 io 

1.36 

2.40 

0.253 

0.459 

4.56 

10.39 

.25 

10.20 

7.80 

1.360 

.496 

1.942 

.61 

• 293 

.383 

. 4 oo 

1.33 

2.52 

.365 

.443 

A. 56 

10.39 

.50 

10.20 

5-40 

1.968 

•599 

2.541 

.66 

.288 

.396 

• 393 

1.25 

2.72 

.565 

.418 

4.56 

10.39 

.75 

10.20 

2.60 

4.088 

.829 

4.607 

.73 

.287 

.417 

• 391 

1-12 

3.34 

1.26 

• 374 

A. 56 

10.39 

.93 

10.20 

.69 

15.293 

.124 

15.252 

.71 

.287 

.431 

• 391 

.95 

6.61 

4.67 

.319 

9-12 

41.58 

.0 

2.55 

41.58 

.255 

.283 

.570 

.91 

.254 

• 295 

• 319 

1.36 

2.20 

.007 

.453 

9-12 

41.58 

.25 

2.55 

31.18 

.341 

.317 

.765 

.93 

.256 

.301 

.312 

1.31 

2.36 

.225 

.436 

9.12 

41.58 

.50 

2.55 

20.80 

•511 

• 359 

1.044 

.70 

.254 

.308 

• 319 

1.26 

2.71 

.502 

.421 

9.12 

41.58 

.75 

2.55 

10. 4 o 

1.022 

.441 

1.602 

.61 

.248 

.320 

.322 

1.18 

3.01 

.959 

.393 

9.12 

41.58 

.93 

2.55 

2.78 

3.823 

.811 

4.352 

.73 

• 249 

• 337 

• 311 

.92 

4.83 

2.97 

• 305 

12.16 

73.92 

.0 

1.43 

7-39 

1.437 

.188 

.315 

.36 

.244 

.262 

.303 

1.48 

.0 

-.078 

.495 

14.59 

106.44 

.0 

1.00 

• 106.44 

•099 

.097 

.228 

’ .27 

.241 

.240 

.298 

1.89 

2.50 

-.123 

.626 

15.20 

115.50 

.25 

.92 

86.62 

.123 

•152 

.277 

.34 

.226 

.246 

• 279 

i.54 

2.46 

.072 

•515 

15.20 

115.50 

.50 

.92 

57.75 

.18b 

.234 

.396 

.38 

• 213 

.260 

.263 

1.32 

1 2.44 

.268 

.438 

15.20 

115.50 

.75 

•92 

28.88 

.368 

.326 

.818 

.87 

.214 

.274 

.264 

1.17 

3.02 

.967 

.391 ■ 

; 15-20 

115.50 

.93 

.92 

7.65 

1.389 

.498 

.195 

.61 

.217 

.286 

.268 

.99 

4.32 

2.44 

.330 

18.24 

166.32 

.50 

.64 

83.19 

.128 

.162 

.286 

.35 

.207 

.237 

.254 

1.44 

2.62 

.301 

.482 

18.24 

166.32 

.75 

.64 

41-59 

.255 

.284 

•570 

.91 

.196 

• 255 

.241 

1.17 

3.03 

.985 

• 391 

18 . 24 

166.32 

.93 

.64 

11.20 

• 948 

.430 

1.527 

.61 

.203 

.269 

.250 

1.00 

4.53 

2.61 

• 334 

21.28 

226.38 

.75 

.47 

56.60 

.188 

.238 

.405 

.38 

.190 

.237 

.234 

1.19 

2.94 

.855 

• 397 

21.28 

226.38 

•93 

.47 

15.18 

• 700 

.392 

1.269 

.63 

.190 

.252 

.234 

.99 

5.02 

3-03 

.330 

24.32 

295.68 

•93 

. 4 o 

19.76 

.538 

.364 

1.079 

.69, 

.181 

•239 

.222 

.97 

5.17 

2.89 

.322 
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TABLE II.- TABULATION OF FLOW CONDITIONS AND PRIMARY RESULTS - Continued 


(l) p e = 1 atm; H w = A. 10 km 2 / sec 2 ; T w = 3000° K; <% = 0 .69 


Uco , 

km/ sec 

H e , 

km 2 / sec 2 

ue 2 /2He 

gwXlO 2 

he j 

km 2 / sec 2 

h-w 

—■ x 10 
He 

le 


a e 


f nnw 

Nu/s/ Hy 

cf \ J&e 

f/s 

y 


A. 56 

10. Ao 

0.0 

39 - Ao 

10.39 

3 - 9^7 

O.6A0 

5.080 

0.61 

0.229 

0.A12 

0.A63 

1.26 

2.76 

0.571 

0.A21 

A. 56 

10. Ao 

.25 

39. Ao 

7.80 

5.262 

.720 

6.157 

.61 

.215 

.A23 

.A32 

1.22 

2.86 

.7A0 

.A07 

A. 56 

10. Ao 

• 50 

39 - Ao 

5.20 

7.891 

.886 

8.295 

.66 

.202 

.A38 

.A08 

1.15 

3-13 

1.06 

.379 

A. 56 

10. Ao 

•75 

39. Ao 

2.60 

15 . 78 A 

.120 

lA.6ll 

•731 

.192 

.A57 

.388 

1.02 

A. 01 

2.01 

.3A0 

A. 56 

10. Ao 

•93 

39. Ao 

.69 

59 -OA 3 

.180 

A8.372 

•71 

.186 

.A73 

.281 

.86 

8.61 

6.76 

.290 

9-12 

Al.58 

.0 

9.70 

A1.58 

.987 

• AlO 

1.809 

•91 

.270 

.3A3 

.367 

1-31 

2.31 

.09A 

.AAl 

9-12 

Al.58 

•25 

9.70 

31.18 

1.316 

.A60 

2.A26 

•93 

.271 

■ 3A9 

.368 

1.26 

2.53 

.331 

.A25 

9-12 

Al.58 

•50 

9-70 

20.80 

1-973 

• 520 

3.313 

• 70 

.268 

•357 

.363 

1.21 

2.76 

.628 

.A08 

9-12 

Al.58 

•75 

9.70 

8.55 

A. 801 

.691 

5.781 

.61 

.258 

.372 

• 350 

I.09 

3 -Ao 

1.3A 

.369 

9.12 

Al.58 

•93 

9.70 

2.78 

1A.761 

-117 

13.803 

• 73 

.256 

.386 

.3A7 

•87 

5.23 

3.37 

.293 

lA .59 

106. AA | 

.0 

3.85 

; 106. AA , 

.385 : 

.1A0 

• 72A , 

• 27 

•273 

.28A' 

.3A8 

1.86 

2.55 

-.076 

.618 

15-20 

! 115.50 
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TABLE II.- TABULATION OF FLOW CONDITIONS AND PRIMARY RESULTS - Continued 
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TABLE II.- TABULATION OF FLOW CONDITIONS AND PRIMARY RESULTS - Continued 
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TABLE II.- TABULATION OF FLOW CONDITIONS AND PRIMARY RESULTS - Concluded 
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Figure 1.- Concluded. 
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Figure 3 .- Continued. 
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Figure 4.- The effect of wall temperature on the correlation parameter K(T W , p e ) 
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Figure 5.- Comparison with other theoretical heat-transfer results. 
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Figure 8.- Comparison of skin- friction parameter with other theoretical 

results . 
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Figure 8.- Concluded. 
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Figure 11.- Boundary-layer displacement thickness; p e = 1 atm, T w - 1000° K. 
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Figure 12.- Boundary-layer momentum thickness; p e = 1 atm, T w - 1000° K. 
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